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Abstract 

The  problem  of  recovering  information  from  single  component  linear  translation  invariant 
systems  is  inherently  ill-posed.  However,  ill-posedness  may  be  circumvented  in  a  multichan¬ 
nel  system  if  the  components  of  the  system  satisfy  the  strongly  coprime  condition.  The 
signal  may  be  completely  recovered  from  a  strongly  coprime  system  by  filtering  the  output 
of  each  channel  with  a  deconvolver,  and  adding.  This  approach  has  been  labeled  multichannel 
deconvolution. 

We  explore  the  state  of  the  art  for  the  multichannel  theory.  We  then  use  sampling  theory 
to  develop  an  alternative  method  for  creating  these  systems.  Modulation  techniques  are 
then  used  to  create  a  strongly  coprime  system,  for  which  the  corresponding  deconvolvers  are 
developed.  We  close  by  discussing  several  applications  of  the  theory. 


1  Introduction 


Linear,  translation  invariant  systems  (e.g.,  sensors,  linear  filters)  are  modeled  by  the  con¬ 
volution  equation  s  =  f  *  p,  where  /  is  the  input  signal,  p  is  the  system  impulse  response 
function  (or,  more  generally,  impulse  response  distribution),  and  s  is  the  output  signal.  We 
refer  to  jj,  as  a  convolver.  In  many  applications,  the  output  s  is  an  inadequate  approximation 
of  /,  which  motivates  solving  the  convolution  equation  for  /,  i.e.,  deconvolving  /  from  p.  If 
the  function  p  is  time-limited  (compactly  supported)  and  non-singular,  we  have  shown  that 
this  deconvolution  problem  is  ill-posed  in  the  sense  of  Hadamard  (see  [16]). 

A  theory  of  solving  such  equations  has  been  developed.  It  circumvents  ill-posedness 
by  using  a  multichannel  system.  If  we  overdetermine  the  signal  /  by  using  a  system  of 
convolution  equations, 


Si  =  f  *  m  ,  *  =  1,...  ,n, 


(1) 


the  problem  of  solving  for  /  is  well-posed  if  the  set  of  convolvers  {/i*}  satisfies  the  condi¬ 
tion  of  being  what  we  call  strongly  coprime.  In  this  case,  there  exist  compactly  supported 
distributions  (deconvolvers) 

Vi,  i  =  1,...  ,n 


such  that 


Hi  ■  V\  + . 

.  .  +  fin  *  Vn  =  1  . 

(2) 

Transforming,  we  get 

Hi  *vi  + . 

•  •  T  Hn  *  vn  =  5 , 

(3) 

which  in  turn  gives 

S\  *  V\  +  . . 

•  •  +  Sn  *  vn  =  /  . 

(4) 

The  theory  of  deconvolution  presented  in  this  paper  is  contained  in  a  larger  group  of 
results  in  the  theory  of  residues  of  analytic  functions  and  their  generalities,  for  example, 
intersection  varieties.  These  results  have  appeared  in  a  series  of  papers  by  Berenstein,  Gay, 
Taylor,  Yger  et  al.  (see  [2]  -  [17]),  and  can  be  interpreted  as  results  in  division  problems, 
interpolation  of  analytic  functions,  analytic  continuation,  digital  to  analog  conversion,  and 
complexity  theory.  For  deconvolution  and  other  applications  to  signal  and  image  processing, 
the  theory  focuses  on  solving  the  analytic  Bezout  equation,  i.e.,  for  given  holomorphic  fi  and 
4>  satisfying  certain  growth  conditions,  solving  for  holomorphic  &  satisfying  the  same  growth 
conditions  such  that 


/l  '  9l  +  ■  ■  •  +  fn  ’  9n  —  4>  •  (5) 

For  our  purposes,  we  want  growth  conditions  given  by  the  Paley-Wiener-Schwartz  Theorem 
and  (j)  —  (f>\,  with  4>\  — »  1  as  A  — >  oo  {<j>\  is  the  transform  of  an  approximate  identity). 
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Solutions  to  Bezout  equations  have  yielded  results  in  deconvolution,  complexity  theory,  so¬ 
lutions  to  systems  of  PDE’s,  theorems  about  interpolation  and  continuation  of  analytic 
functions,  and  results  in  number  theory  (see  [2]  -  [17]). 

We  describe  in  section  2.1  the  strongly  coprime  condition,  and  we  give  examples  of  sets 
of  strongly  coprime  system  response  functions  and  their  deconvolutions  for  functions  in  one 
and  several  variables.  In  the  language  of  applications,  the  set  of  convolvers  {/ij}  models  a 
linear  translation  invariant  multichannel  system  consisting  of  an  array  of  sensors  or  filters. 
The  system  is  created  so  that  no  information  contained  in  the  input  signal  /  is  lost.  The 
signal  /  is  gathered  by  this  system  as  {s*  =  /  *  //*}.  The  signals  s*  are  then  filtered  by  the 
Vi  (which  have  been  created  digitally,  optically,  etc.,  in  coordination  with  the  creation  of  the 
system  and  possibly  tailored  to  be  optimized  under  some  constraint)  and  added,  resulting  in 
the  reconstruction  of  /.  We  discuss  the  various  classes  of  impulse  response  functions  modeled 
by  the  theory. 

Section  2.2  gives  the  development  of  a  real-variable  method  of  solving  Bezout  in  a  specific 
case.  Using  Shannon  sampling  in  the  frequency  domain,  this  development  utilizes  the  zero 
sets  of  the  {&}  as  two  different  sampling  rates.  Section  3  contains  a  new  method  for 
creating  strongly  coprime  systems.  In  particular,  a  given  fixed  convolver  is  modulated  a 
certain  amount,  resulting  in  a  strongly  coprime  pair.  The  deconvolvers  associated  with  this 
system  are  then  created.  This  approach  gives  a  new  technique  for  constructing  multichannel 
systems.  Section  3.2  gives  a  development  of  these  deconvolvers  using  sampling. 

We  close  in  section  4  with  a  discussion  of  applications  of  the  multichannel  theory  to 
signal  and  image  processing.  In  a  system  designed  to  do  signal  or  image  processing,  the 
deconvolvers  could  act  as  preliminary  image  or  signal  enhancement  operators,  restoring  high 
frequency  information.  These  operators  can  work  simultaneously  with  other  operators  which 
analyze  the  image  or  signal,  for  example,  wavelet  operators  (see  [16]).  The  work  with  wavelet 
analysis  in  image  processing  looks  particularly  promising.  Deconvolution  will  restore  the  high 
frequency,  fine-detail  of  the  image.  Wavelet  methods  can  then  be  used  for  image  compression 
and  analysis,  e.g.,  for  edge  detection,  motion  detection,  feature  extraction,  etc. 

We  need  some  background  information.  For  an  integrable  function  /,  we  define  the 
Fourier  transform  f(u)  as 

/M  =  [  f(t)e~2^dt, 

JRn 

for  t  £  W1  (time),  u>  £  Rn  (frequency)1.  Given  an  integrable  function  g  on  Rrt,  the  inverse 
transform  is 

g(t )  =  L  g(u))  e27rt^  dio  . 

J  R" 

If  f(t)  is  a  compactly  supported  function,  f(u)  is  an  analytic  function.  Moreover,  f(u>) 
analytically  continues  to  all  of  n-dimensional  complex  Euclidean  space  O1.  This  continuation 
gives  the  Fourier-Laplace  transform  of  /,  which  is  defined  by 

f(0  =/  f(t)e~2^dt,  (GC. 

JRn 

1  We  define  the  Fourier  transform  in  this  way  so  that  ||/||2  =  ||/||2-  Some  “standard  formulae”,  e.g., 
Shannon  sampling  as  presented  in  [28] ,  have  been  appropriately  adjusted. 
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(See  [19],  [24],  [28],  [29]  for  further  reference.) 

The  process  of  deconvolution  uses  some  of  the  basic  theory  of  distributions.  Let  N  denote 
the  natural  numbers,  and  let  Z  denote  the  integers.  We  define  three  sets  of  test  functions, 

V(Rn)  =  C^°(Rn)  =  {(j)  e  C°°(Rn)  with  compact  support} , 

dk 

<S(Rn)  =  {(j)  G  C0C(Rn)  with  lim  tn—rMt)  —  0}  ,  for  all  n,  k  e  N,  and 

|t|— kx>  atK 

£(Rn)  =  C°°(RL). 

Each  of  these  has  an  associated  space  of  distributions,  denoted  £'(Rn),  <S'(Rn),  and  X>'(Rn). 
The  action  of  a  distribution  T  on  a  test  function  (j>  is  denoted  by  (T;  <j>).  The  more  restrictive 
the  notion  of  convergence  on  the  space  of  test  functions,  the  broader  the  class  of  distributions. 
We  have  that 

£'  C  S'  C  V' . 

The  Fourier-Laplace  transform  is  defined  on  S'  and  £',  the  classes  of  tempered  and  compactly 
supported  distributions,  respectively.  The  sets  of  transforms  are  denoted  by  S'  and  £'.  For 
A  €  £',  the  Fourier  transform  continues  as  an  analytic  function  to  all  of  C1 .  This  continuation 
defines  the  Fourier-Laplace  transform  on  £',  which  is  given  by 

A  =  ^A;  e~27"^>) . 

Properties  of  the  function  A  are  given  in  the  following. 

Theorem  1.1  (Paley-Wiener-Schwartz  (PW5)[1])  a.)  The  Fourier-Laplace  transform 
of  an  infinitely  differentiable  function  f  with  compact  support  C  (|/.|  <  A}  C  Rn  is  an  entire 
function  f(()  in  Cn  which  satisfies  the  following  property: 

(Pi)  For  every  integer  N  >  0  we  can  find  a  positive  constant  C  —  C(N)  such  that 

|/(C)|  <  C(  1  +  KirV^^I  for  all  (eC.  (6) 

Conversely,  every  entire  function  in  CL  satisfying  property  (Pi)  is  the  Fourier-Laplace 
transform  of  a  C°°  function  with  compact  support  C  {|t|  <  A}  in  Rn. 

b.)  The  Fourier-Laplace  transform  of  a  distribution  A  with  compact  support  C  {\t\  < 
A}  C  Rn  is  an  entire  function  X(()  in  CL  which  satisfies  the  following  property: 

(P2)  There  is  a  positive  constant  C  and  an  integer  N  >  0  such  that 

|a(C)|  <  ^(1  +  iCD^e2^1^1  for  all  C  e  CL.  (7) 

Conversely,  every  entire  function  in  CL  satisfying  property  (P2)  is  the  Fourier-Laplace 
transform  of  a  distribution  A  with  compact  support  C  {\t\  <  A}  in  Rn. 

(See  [1],  [24]  for  further  reference.) 
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2  Deconvolution  from  a  System  of  Convolution  Equa¬ 
tions 

2.1  An  Overview  of  the  Theory 

We  begin  by  giving  a  framework  of  realistic  deconvolution  problems.  For  these  problems, 
our  input  data  is  p  and  /  *  p,  and  our  proposed  solution  is  /.  For  a  fixed  p,  we  have 
an  associated  convolution  operator  C^(/)  =  /  *  p .  If  is  injective,  then  the  inverse,  or 
deconvolution  operator,  is  *  p)  =  / .  The  deconvolution  problems  we  consider  are 

for  convolvers  p  which  are  realistic  mathematical  models  of  the  impulse  response  functions 
of  linear  translation  invariant  systems.  Therefore,  we  exclude  distributions  which  are  not 
compactly  supported  (since  one  would  have  to  integrate  for  all  time  to  get  any  information 
from  such  a  system),  distributions  of  order  k  >  1  (since  one  would  have  to  impose  smoothness 
conditions  on  any  input  functions),  and  any  measures  that  are  singular  with  respect  to 
Lebesgue  measure  (since  such  a  system  is  impossible  to  build).  If  p  is  realizable,  then  the 
mapping  C Lp{ Rn)  — >  1 Z^  C  Z/(Rn),  where  1Z^  is  the  range  of  is  continuous  and 
linear.  The  most  general  convolver  we  will  consider  is  a  “well-behaved”  measure  as  defined 
below. 

Definition  2.1  The  distribution  p  is  a  realizable  convolver  if  p  is  a  compactly  supported 
finite  Borel  measure  which  is  absolutely  continuous  with  respect  to  Lebesgue  measure  on  Rn . 

The  Radon-Nikodym  Theorem  gives  us  that  all  realizable  deconvolvers  act  as  compactly 
supported  Ll  functions.  The  set  of  realizable  convolvers  also  includes  the  set  of  non-singular 
probability  measures  of  compact  support  n.  Given  such  a  measure,  the  Radon-Nikodym 
derivative  is  the  probability  density  function  of  7T. 

We  have  shown  the  following. 

Theorem  2.1  ([16])  Let  p(t)  be  a  realizable  convolver.  Then  for  f  G  C(Rn),  the  con¬ 
volution  operator  CM(/)  =  /  *  p  is  not  injective.  Therefore ,  the  deconvolution  problem  of 
recovering  f  from  f  *  p  is  ill-posed  in  the  sense  of  Hadamard. 

Remark  :  The  method  of  proof  for  the  theorem  above  gives  an  indication  of  the  “degree  of 
ill-posedness”  for  inverting  a  given  convolution  operator.  The  procedure  involves  examining 
the  set  of  all  linear  combinations  of  functions  g  such  that  supp(<?)  C  Z^,  where  Z^  is  the 
variety  defined  by  the  zero  set  of  fi,  the  Fourier-Laplace  transform  of  p.  If  p  is  realizable,  Z^ 
is  not  empty.  Functions  g  as  above  will  be  in  the  kernel  of  the  operator.  The  signals  modeled 
by  these  functions  will  all  be  undetectable  by  the  sensor  modeled  by  p.  For  example,  in  one 
variable,  given  any  Co  =  +  w? o  €  Z —  C0  G  Z^,  and  so  g(t)  —  e-2^0*  cos(2-7r co0t)  ,  h(t)  = 
e— 2ir»?ot  sjn (2nco0t)  are  functions  such  that  supp(<?(C))  and  supp(h(C))  are  contained  in  Z^. 
Moreover,  as  there  are  infinitely  many  pairs  (C,  — C}  of  zeros  in  Z/x>  there  are  infinitely  many 
such  functions  (see  [16]). 

If  /  G  L2(Rn),  we  have  that  L2(Rn)  — >  IZ^  C  L2(Rn)  is  injective  with  ||/  *  p\\i  < 
ll/ltalMI-  We  have  shown  that  the  deconvolution  operator  IZ^  — >  L2( Rn),  defined  by 
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X>M(/*/x)  =  /,  is  discontinuous.  We  do  this  by  building  a  sequence  {/„}  of  L2(Rn)  functions 
such  that  1 1 /n 1 1 2  — >  oo  as  n  — >  oo,  but  |j/n  *  p H2  — »■  0  as  n  — >  00,  for  any  fixed  realizable 
convolver  p.  Thus,  this  deconvolution  problem  violates  the  third  condition  of  well-posedness. 

Theorem  2.2  ([16])  Let  p  be  a  realizable  convolver.  Then  for  f  G  L2(Rn),  the  deconvo¬ 
lution  operator  V^f  *  p)  =  /  is  an  unbounded  and  therefore  discontinuous  linear  opera¬ 
tor.  Thus,  the  deconvolution  problem  of  recovering  f  from  f  *  p  is  ill-posed  in  the  sense  of 
Hadamard. 

These  theorems  generalize  in  the  following  way.  If  p  is  realizable,  and  /  €  Lp( Rn),  then 
the  deconvolution  operator  T>n(f*p)  —  f  is  an  unbounded  and  therefore  discontinuous  linear 
operator. 

We  have  constructed  solutions  to  the  deconvolution  problem  for  various  classes  of  com¬ 
pactly  supported  convolvers.  We  have  assumed  only  that  our  input  or  initial  functions  / 
are  of  finite  energy.  Ill-posedness  is  circumvented  by  creating  a  multichannel  system.  Each 
channel  in  the  system  is  a  convolver,  and  the  system  overdetermines  the  input  function  with 
all  of  the  convolvers  in  an  array  chosen  so  that  any  information  lost  by  one  convolver  is 
retained  by  another.  The  theory  of  deconvolution  discussed  in  the  following  paragraphs  has 
its  roots  in  the  work  of  Wiener  [35]  and  Hormander  [25],  and  has  been  developed  into  a 
working  theory  by  Berenstein,  Taylor,  Gay,  Yger,  et  al.  [2]  -  [17].  These  methods  are  both 
linear  (convolution  with  deconvolvers)  and  realizable  (the  support  of  the  deconvolvers  being 
contained  in  the  bounded  support  of  the  kernels  of  the  convolution  equations).  Thus,  de- 
convolution  at  a  point  t  €  Rn  depends  only  on  data  near  t.  The  theory  assumes  no  a  priori 
information  about  the  input  signals.  Moreover,  the  theory  can  be  used  to  develop  a  stable 
system  for  complete  signal  recovery. 

The  theory  starts  with  the  following  existence  theorem  of  Hormander. 

Theorem  2.3  (Hormander  [25])  For  the  compactly  supported  distributions  {//,}£_  x  on 
Rn,  there  exist  compactly  supported  distributions  1  such  that 

pi  *  i/i  +  ...  +  pn*  vn  =  6 

if  and  only  if  there  exist  positive  constants  A  and  B  and  a  positive  integer  N  such  that 

(£  IRK)I2)  ’  >  +  ICir" .  C  e  C* .  (8) 

Definition  2.2  A  set  of  convolvers  {/^}”=1  that  satisfy  the  inequality  in  the  theorem  is  said 

to  be  strongly  coprime. 

We  note  that  the  only  way  a  single  compactly  supported  convolver  can  be  strongly 
coprime  is  for  it  to  be  a  translate  of  the  identity  convolver;  that  is,  for  it  to  be  the  Dirac 
delta  or  a  translation  thereof.  This  result  is  one  way  to  state  the  general  ill-posedness  of 
single  convolution  equations  under  the  constraints  imposed  by  the  conditions  of  the  theorem. 
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Definition  2.3  Let  {Fi}f=1  be  a  given  set  of  functions  in  £'(Rn).  A  solution  to  the  analytic 
Bezout  equation 


5><(C)G4(0  =  1  (9) 

i 

is  a  set  {(rj}f=1  in  £'(Rn)  that  satisfies  the  equation. 

By  PWS  and  basic  properties  of  the  Fourier- Laplace  transform,  a  solution  to  the  analytic 
Bezout  equation  is  equivalent  to  solving  for  a  set  {t/J  in  £'  such  that  YiVi  *  »i  =  5,  for  a 
given  set  {/z,}  in  £'.  By  Hormander’s  theorem,  a  strongly  coprime  set  {//,}  is  precisely  a 
set  for  which  the  analytic  Bezout  equation  has  a  solution.  The  strongly  coprime  condition 
guarantees  not  only  that  the  transforms  of  the  convolvers  have  no  common  zeros,  but  also 
that  these  zeros  do  not  cluster  too  quickly  as  |£|  — >  oo.  Thus,  if  a  given  signal  /  is 
overdetermined  by  a  strongly  coprime  system  of  convolution  equations, 

Si  =  f  *  Hi  ,  i  =  1, . . .  ,  n , 

then  the  problem  of  solving  for  /  is  well-posed.  We  solve  for  a  set  {^(C)}  of  deconvolvers 
which  satisfy  the  analytic  Bezout  equation 

YtiWi  =  1- 

i 

Taking  inverse  transforms  of  both  sides  of  Bezout  gives 

Y, Vi  *"i  =  5- 

i 

The  fact  that  the  strongly  coprime  condition  is  the  inversion  of  the  PWS  growth  bound 
allows  us  to  solve  for  deconvolvers  that  are  compactly  supported.  This  in  turn  yields 


Si  *  V\  +  . . .  +  sn  *  un 
=  (/  *  Ail)  *  Vi  +  ...  +  {f  *Vn)*  Vn 

=  f  *  (vi  *  J'l)  +  •  •  •  +  /  *  (Vn  *  Vn) 

=  f  *  (Vl  *  "1  +  •  •  •  +  Vn  *  Vn) 

=  f*5 
=  /• 

Thus,  the  deconvolution  problem  can  be  solved  by  constructing  the  Dirac  6  for  a  given  class 
of  convolvers.  This  construction  begins  with  a  solution  to  the  analytic  Bezout  equation.  The 
inverse  transforms  of  these  solutions  are  compactly  supported  distributions. 

To  construct  compactly  supported  deconvolving  functions,  we  begin  by  solving  a  more 
general  analytic  Bezout  equation,  i.e.,  for  given  analytic  fi]  and  '0  satisfying  PWS  growth 
conditions,  solving  for  analytic  9i  satisfying  PWS  growth  conditions  such  that 

fi[  ■  v{  +  .. .  +fiPn  •££  = 
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Moreover,  we  want  0  =  0a  >  with  0a  — >  1  as  A  — >  oo  (0a  is  the  transform  of  an  approximate 
identity).  This  gives  us  deconvolving  functions,  i.e.,  deconvolvers  {z^}  such  that  Hi*t/i,^>  + 
...+  fin*  vn^  =  0,  which  in  turn  give 


(/  *  fli)  *  Vij,  +  .  .  .  +  (/  *  fin)  *  Vn,*  =  f  *  0  =  U  ■ 


Then,  as  0  — ¥  5,  — ¥  /  in  the  sense  of  distributions.  The  deconvolvers  in  these  imple- 

mentable  formulae  are  periodic  functions  expressed  in  their  Fourier  series  expansions. 

We  give  an  example.  Let  tel,  p  be  a  prime  number,  and  let 


Hl(t)  —  X[_l)l](i)  ,  —  '^■[~y/p,y/p\{t) 


model  the  impulse  response  of  the  channels  of  a  two-channel  system.  Thus 

- — -  /  ^\  sin(27rC)  sin(27rv^O 

=  ■ 


Let 


k  e  N 


denote  the  zero  sets  of  /Ti(C),  /£?(())  respectively.  An  examination  of  the  Fourier-Laplace 
transforms  /£((),  i  =  1,2,  gives  that  {/a}  is  strongly  coprime  (see  [6]).  We  choose  an 
arbitrarily  close  approximation  0  of  the  Dirac  6  based  on  certain  criteria,  i.e.,  0  in  C 4  with 
support  in  (—(1  +  y/p),  (1  +  y/p))-  Then  |0(z)|  <  for  z  £  Zi  U  Z2.  The  smoothness 

and  the  size  of  the  support  of  0  guarantee  that  the  deconvolvers  are  compactly  supported. 


Theorem  2.4  ([16])  The  set 


pi(t)  =  X[_i,i](i)  ,  p2(t)  =  X[_^>V5](t) 

is  a  strongly  coprime  pair  of  convolvers.  Let  f  G  L2(K)  and  let  0  be  a  C4  function  with 
support  in  (—(1  +  y/p),  (1  +  y/p))  such  that  0  >  0  and  /f^,0(t)  dt  =  1.  The  deconvolvers 
Vi ^  such  that 

f  *  0  =  (/  *  hi)  *  vl>xjj  +  (f*  H2)  *  ^ 

are  given  by  the  formulae 


vi  At)  = 


**,*(*) 


Y'  ( l  (el™(t+y/p)  _ 

zTz,  Mi (z)ftH2(z)  ^  V 

(]_  (e2mz(t+l)  _ 

fk,  fc{z)icHi{z)  \z  ' 


(10) 


(11) 


The  function  /  *  0  is  an  arbitrarily  close  approximation  of  f  which  converges  to  f  in  the 
sense  of  distributions  as  supp(0)  — ¥  {0}. 
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The  deconvolvers  in  Theorem  2.4  were  developed  by  using  the  Jacobi  interpolation  for¬ 
mula  and  the  Cauchy  residue  calculus  in  the  complex  plane  (see  [16]).  It  is  also  possible 
to  use  real- variable  methods,  in  particular  Shannon  sampling,  to  create  the  This  is 

done  in  section  2.2. 

This  development  works  for  other  classes  of  convolvers  and  filters.  The  current  stock 
of  convolvers  and  their  associated  deconvolvers  includes  characteristic  functions  of  squares 
and  (hyper)cubes  (see  [3],  [13],  [16]),  and  characteristic  functions  of  disks  and  n-dimensional 
balls  (see  [4],  [5]).  The  theory  has  been  expanded  in  one  variable  to  more  general  convolvers, 
including  convolvers  modeled  by  linear  combinations  of  characteristic  functions,  linear  com¬ 
binations  of  n-fold  convolutions  of  characteristic  functions  with  equally  space  knots  (cardinal 
splines),  and  truncated  sine,  cosine,  and  Gaussian  functions  (see  [17]).  We  have  shown  the 
conditions  for  a  strongly  coprime  set  of  convolvers  {//*}  for  each  of  these  types  of  functions, 
which,  as  in  the  example  above,  is  a  condition  on  the  zero  sets  of  the  Fourier-Laplace  trans¬ 
forms  {Jm}.  We  then  have  solved  for  deconvolvers  {u^}  such  that  pi*v\^-\-. . 
where  if  is  an  approximate  convolution  identity,  by  solving  the  modified  Bezout  equation 
Ah  ’  ^1  ,i/>  T  •  .  .  T  Hn  '  =  0- 

In  several  variables,  the  formulae  for  the  deconvolvers  is  simplified  by  not  only  solving  for 
an  approximation  to  the  6,  but  also  by  strengthening  the  strongly  coprime  condition.  Let 


2  =  {z  G  O1  :  px(z)  =  . . .  =  Hm(z)  =  0}  . 


We  say  that  Z  is  almost  real  if  there  exits  a  constant  a  >  0  such  that 

2  C  [z  g  r  :  |S*|  <  alog(2  +  |*|)}  . 

An  almost  real  set  is  discrete  (see  [3]).  Let  Br  =  {|z|  <  r}.  Define  the  counting  function 

n(Z,  r )  =  card  {2  n  Br}  . 


Also  define  the  distance  function 

d(z,  2)  =  min(l,  min  {\z  -  C|  :  C  G  Z}) . 

Let  Hi(9)  denote  the  supporting  function  of  the  convex  hull  of  Usupp^i,  i.e.,  for  9  G  Rn, 

Hx  =  maxj  max  {(x,  9)  :  x  G  supp  fij}  . 

Definition  2.4  ([3])  A  family  ofm  distributions  of  compact  support  in  R”  is  well  behaved 
if  there  exist  positive  constants  A,  B,  N ,  k  and  a  supporting  function  Hq  with  0  <  Hq  <  Hi 
such  that  2  is  almost  real, 

n(Z,r)  =  0(rA) , 

and 

(m  \  2 

E  IS(C)|!)  >  Bd(z,  2)’eH“<32)(  1  +  |(|)- "  .  (12) 
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Definition  2.5  ([3])  A  well  behaved  system  {m}iL  1  is  called  very  well  behaved  if  there 
exist  positive  constants  f3  and  M  such  that  for  every  z  e  Z,  we  have 

\J(z)\>  P(l  +  \z\)~M ,  (13) 


where  \J(z)\  = 


det 


djiA 
dzj . 


For  a  very  well  behaved  system,  the  zeros  of  the  pLi  are  simple.  For  convolvers  in  n- 
dimensional  Euclidean  space  Rrl ,  we  have  that  a  very  well  behaved  system  needs  only  n  +  1 
elements. 


Theorem  2.5  ([3])  Given  a  very  well  behaved  system  {/Xj}^1  and  ip  «  C°°  function  with 
compact  support  in  Rn ,  there  exist  compactly  supported  functions  such  that  for  all 

(eC 


II 

#(*) 

"}((,z)  .. 

.  »r+i(c,*) 

J(z)fXfXi{z) 

Ph(C>z)  .. 

Mi(0  ■ 

■  Pn+l{C) 

(14) 


where 


H(C,Z)  = 


1^3  (Cl)  •  •  ■  ;  Cki  %k+ 1)  •  •  •  i  Zn)  Mj(Cl)  •  •  •  >  Cfc— X;  Zki  •  •  •  i  Z-n) 

(Cfe  -  Zk) 


The  function  ip  can  be  an  arbitrarily  close  approximation  of  the  6  which  converges  to  the  5 
in  the  sense  of  distributions  as  supp(-0)  — >  {0}. 


The  formula  given  in  the  theorem  above  gives  a  solution  to  a  generalization  of  the  Bezout 
equation  for  a  very  well  behaved  family.  The  function  ip  can  be  chosen  arbitrarily  close  to 
the  6.  Again,  the  formulae  for  the  deconvolvers  are  functions  expanded  in  Fourier  series,  not 
distributions. 

This  formula  gives  the  solution  to  the  deconvolution  problem  when  the  convolvers  are 
(hyper) rectangular  regions  in  Rn.  It  has  a  variety  of  other  applications,  including  producing 
a  more  explicit  version  of  the  Ehrenpreis  Fundamental  Principle,  i.e.,  the  representation  of  all 
the  solutions  of  a  system  of  linear  partial  differential  equations  with  constant  coefficients  [12], 
characterizations  of  interpolation  varieties,  and  new  versions  of  gap  theorems  for  Dirichlet 
series  [12]. 

Certain  problems  on  deconvolution  can  be  phrased  as  Pompeiu  problems.  We  give  the 
following  example.  Let  W,  X2  be  the  characteristic  functions  of  the  disks  S(0,  ri),  £(0,r2), 
and  let  E  be  the  collection  of  positive  quotients  of  zeros  of  the  Bessel  function  J\. 


Theorem  2.6  ([4])  Let  f  be  a  continuous  function  in  the  plane.  If  there  exists  an  a  >  0 
such  that  >  \rj~a  for  all  £,r)  e  E  with  £,77  >  0,  then  the  mapping  P:  f  — » 

(%!  *  /,  X2  *  /)  is  injective.  Moreover ,  there  exist  iq,  u2  such  that 

v\  *  (W  *  /)  +  ^  *  (X2  *  /)  =  / .  (15) 
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In  [4],  this  theorem  was  also  extended  to  the  local  problem  of  reconstructing  /  in  some 
disk  B(0,  R),  R>  r\+  r2,  from  its  averages  on  B( 0,  ri),  B(0,  r2). 

The  result  fits  in  the  larger  context  of  Pompeiu  problems.  Let  Ex, . . .  ,  Em  be  compact 
sets  of  positive  measure  in  R",  let  C(Rn)  denote  the  space  of  continuous  functions,  and 
let  M(n )  be  the  group  of  Euclidean  motions  in  R”.  Then,  the  (global)  Pompeiu  trans¬ 
form  associated  to  the  sets  Pi,...  ,En  is  the  mapping  P:  P(Rn)  — >■  C(M(n))m  given  by 

(Pf)(t)  =  (Ufdx . f,Em  f  dx)  .  The  problem  is  then  to  give  conditions  on  the  sets 

E\, . . .  ,  Em  to  guarantee  that  the  mapping  P  is  injective.  The  problem  also  has  a  local 
formulation  (see  [5]).  We  can  construct  the  inverse  mapping  by  constructing  deconvolvers 
which  recover  the  function  /,  as  in  the  theorem  above. 

The  application  of  the  theory  from  Pompeiu  problems  to  deconvolution  is  powerful.  These 
techniques  will  allow  us  to  create  strongly  coprime  systems  by  varying  a  single  element.  We 
explore  this  idea  further  in  Section  3. 


2.2  Multichannel  Deconvolution  and  Sampling 

We  now  give  a  construction  of  the  deconvolvers  in  Theorem  2.4  using  the  Shannon  Sampling 
Theorem  in  the  frequency  domain.  First  recall  the  Shannon  Sampling  Theorem. 


Theorem  2.7  (Shannon  [31])  Let  f  be  a  function  of  finite  energy  onR(f  E  L2(R))  with 
Fourier  transform  f(co)  =  0  for  all  |w|  >  LI,  i.e.,  f(t)  is  Ll-band-limited. 


a.)  J/T  <  1/20,  then  for  all  t  E  R, 


OO 

f(t)=T  £  /(n T) 

Tl  —  —  OO 


sin (Jf(t  -  nT )) 
n(t  —  nT) 


(16) 


b.)  7/T  <  1/20  and  f(nT )  =  0  for  all  n  E  Z,  then  f  =  0. 

Functions  in  L2{ R)  that  are  O-band-limited  are  the  transforms  of  analytic  functions  of 
finite  energy  which  satisfy  Paley-Wiener  growth  estimates.  Such  functions  are  said  to  belong 
to  the  set  PWq(R).  The  Shannon  Sampling  Theorem  has  a  dual,  in  that  a  given  O-time- 
limited  function  can  be  sampled  in  frequency  and  reconstructed  (see  [28]). 

To  give  our  alternate  formulae  for  the  approximate  deconvolvers  v^,  we  use  the 
Jacobi  interpolation  formula  to  get  a  solution  to  the  modified  Bezout  equation.  We  then  show 
that  these  deconvolvers  satisfy  certain  PW  growth  estimates.  This  in  turn  allows  us  to  apply 
sampling.  We  use  the  zero  sets  of  the  transforms  of  the  convolvers  to  give  us  the  sampling 
rates.  More  general  modifications  of  this  construction  do  not  use  Jacobi  interpolation  (see 
[32,  33]). 

Theorem  2.8  Let  fxx (t)  =  X[_1;1](t),  jx2(t)  =  for  p  prime.  Let  f  E  L2(R)  and 

let  -0  be  an  even  CA  function  with  support  in  (—(1  +  y/p),  (1  +  y/p))  such  that  0  >  0  and 
ff°0 0  ip(t)  dt  =  1.  The  deconvolvers  such  that 

/  *  0  =  (/  *  jUi)  *  yx^  +  (/  *  Hi)  *  ^ 
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are  also  given  by  the  formulae 


„  m  _  I  _L  V(_1  v+1  ,  _J_  v  *  ,»[./#  I  v 

l'*()  ',2VP,5(  ’  a50/2v®  2y/pHoMn/2y/p)  '  ’ 


(17) 


j.  /UN  |  1  V^/  ili+1  ^OV2)  |  ^  ritint  I  y 

^W-|5E(-1)>  ^72y+2£^V2)e  lXH,H 


(18) 


Proof  :  Let  w  €  R.  Since  0  is  even,  0(a))  is  real.  Consider  first  iq^.  The  formula  for  v2^ 
will  follow  similarly.  Recall  that 


,  k  e  N 


denote  the  zero  sets  of  //2(C))  respectively.  By  equation  (1)  from  Theorem  3.3  in  [16], 

we  have 


"w  =  E 


0(*)  /WO 


(*). 


^(z)tc^(z)  VC  - 

(A  similar  development  is  given  in  Section  3  -  see  equations  (21),  (22).)  We  claim  that 
V\ ^  G  PW^(R).  To  see  this,  first  note  that 


£ 

z£Z2 


0(a) 


Hi(z)-^Li2(z) 


<  00 . 


Now,  for  each  z  G  Z2  and  |C  —  z\  >  T)  >  0  with  97  <  1/2 y/p, 

\M0 


\(~z 

since  ) uft  G  PW’v^(l).  If  |£  —  z\  <  77,  then 

\MO 


<  -|/25(C)I  <  -e27r^|scl 

V  V 


<  sup 


K~z 

Therefore,  there  exists  C  >  0  such  that 


Mz(C) 


C-z 


<  Ce27rv^^^ 


independently  of  z  G  Z2.  Thus, 


V  2(C) 

ST 

$(z) 

1  —  sup 

zez-2 

C-Z 

2s 

zez2 

MZ)^MZ) 
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Note  that  for  each  z  e  Z2, 


[  <  / 

Jr  C  —  z  J\c~. 


C-z 


\M0 


\C-z\>V  C  —  Z 


di + L 


MO 


C-z 


d( 


<  +  n\\±M 


2 

oo 


which  is  bounded  independently  of  z.  Therefore, 


/  \ViAO\2dC<  sup  / 

«/R.  z £ Z'2  “  ® 


MO 


C  —  z 


dc  E 


k2£Z2 


$(z) 


Mz)fcMO 

Now,  since  Z2  =  j|y=}  ,  k  €  N,  and  since  for  n,  me  Z  \  {0}, 


<  oo . 


Mm/2  Vp) 


(m/2^/p)  -  (n/2^/p) 


=  6(n  —  m ) , 


we  have  that  for  n  e  Z  \  {0}, 


-  (JL.\  -  0(n/2VP) 

Ul^2y/p’  JTi{n/  2y/p) 


and 


M°)  =  E 


_ *00 '  !2y/P) _ ^2(0)  _  y^/  jW+l  0(^/2V^) 

^oMj/2Vp)dcMj/2Vp)^~n/2Vp)  m  MJ/2Vp)  ' 


Therefore,  by  Shannon’s  sampling  formula  we  may  write 

=  f  E( 


1^+i  i>U/2Vp)  ^  sin^Try^Q  +  ^  ^n+i  $(n/2y/P)  sin7r(2v^C  ~  ») 


Wo  Mj/2Vp)  )  ^VpC  '  ^  ^  MW2Vp)  7r(2 VpC~u) 

Taking  inverse  Fourier  transforms  and  simplifying  gives 


,  ,,,  (  1  V(  ,u^{i/2VP)  1  * (n/2^) 


Letting  p  =  1  in  the  above  arguments  gives 


^(t) 


2  2 
z  j¥0 


$(n/2) 

/^(i/2)  1  2^0p5(n/2) 


,7r  int 


X[_1;1].  □ 


Figures  1  and  2  [20]  give  simulations  of  this  result. 
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Figure  1.  C°°  Function  (T)  Approximate  Deconvolver  System. 

l.a.)  Channel  1  Deconvolver  V\. 

l.b.)  Channel  2  Deconvolver  v2. 

l.c.)  *  n i  *  v\. 

l.d.)  'I '  *  v2. 

l.e.)  ’T  =  ^  *  Hi  *  Vi  +  tE'  *  fji2  *  v2,  the  resulting  approximate  delta, 
l.f.)  The  Fourier  transform,  or  frequency  response,  of  this  system. 


Figure  2.  Transfer  Response  for  the  System. 

2. a.)  The  response  of  each  channel,  in  the  two  channel  system  of  Figure  1,  to  two  adjacent 
Gaussian  pulses  (denoted  as  input  signal  /). 

(/•h  *  /  —  solid  ,  1^2*  f  ~  dotted) 


2.b.)  Transfer  response  for  the  system  given  in  Figure  1.  Here  the  dotted  line  represents  the 
input  function  while  the  solid  line  shows  the  system  response. 


summ  150 


1e.  =  T*//,  *  v,  +  VF*^2  *  v2 ,  the  resu 


log  ( |  specn 


If.  4/  =  'F 


2a.  The  response  of  each  channel,  in  the  two  channel  system  of  Fig.  1 ,  to  two  Gaussian 
pulses  (denoted  as  the  input  signal,  f). 


2b.  System  transfer  response,  for  the  two  channel  system  of  Fig.  1 ,  when  the  input  function  is  the 
sum  of  two  time  shifted  Gaussian  pulses.  The  dotted  line  is  the  input  function,  the  solid  line  is  the 
deconvolved  function. 


Fig  2.  Transfer  response  for  the  C"  system  shown  in  Figure  1 . 


To  reconcile  the  apparent  differences  between  (10),  (11)  and  (17),  (18),  we  compute  as 
follows.  Let  T\  =  1  and  r2  =  ^/p.  Then 


Thus, 


Hi(n/2rj) 


2sin(n7r(rj/rj)) 

(nn/rj) 


d 

d( 


Hj{n/2rj) 


2  r,'  cos(n7r) 

(n/2  rj) 


(n/^rj) 


(t) 


V"'  fe)  fe2irw(t+r,)  _  -A 

i?,  £(*)*$(*)  V*  ^  J 

y-  fin/2rj)  (n/2rj)  (  1  / 

^0/^(n/2rj)2rjcos(n7r)  V(«/2L>)  ' 


g2tri(n/2r j)(t+rj)  _ 


1  y-,  •0(n/2rJ) 
2ri  ffinfa) 

1  V)(n/2 rj) 

2rj  Mn/2rj) 


(— l)n  ((eiV^')‘  -  l)  *[_,„,,](*) 

((_!)«+! +e«(n/r,)t)  X[_r.)r.]W< 


The  formula  for  Vj^  is  given  by  a  permutation  of  the  indices. 
We  can  simplify  these  formulae  by  noting  that 

j&i(0)j?T(0)  +  JT2( 0)2$(0)  =  -0(0) . 


(19) 


Since 


limpiCC)  =  2  u, 


(19)  will  also  hold  if 


This  gives 


ui,i>  (t) 


2^(0)  = 


^(0) 

4rj 


M  +  y  j(n/2ri) 

4r<  ^0  »i(n/2ri) 


Again,  the  formula  for  is  given  by  a  permutation  of  the  indices.  However,  if  Pj(0)  = 
the  deconvolvers  are  not  continuous  at  ir*.  There  is  a  subtlety  in  this  development.  The 
two  lattices  2X  and  Z2  overlap  at  z  =  0.  To  uniquely  determine  the  we  must  have  either 
continuity  at  the  endpoints  or  must  specify  z/^(0)  and  (see  [33]).  The  deconvolvers 

given  by  (17)  and  (18)  are  continuous. 
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3  Creating  a  System  Using  Modulation 

We  can  use  modulation  to  create  a  strongly  coprime  system.  This  new  technique  for  creating 
these  systems  allows  for  a  greater  flexibility  in  the  development  of  actual  systems.  The  system 
is  created  by  making  two  identical  copies  of  a  given  sensor,  splitting  the  signal  into  two  the 
two  separate  channels,  and  appropriate  modulation  of  both  the  input  and  output  of  one 
of  the  channels.  These  two  outputs  are  then  convolved  with  the  appropriate  deconvolving 
filters  and  added,  resulting  in  the  reconstruction  of  the  input  signal. 

Let  Hi  model  the  impulse  response  of  a  given  system.  Let  E^0  =  e 27TK°L  We  will  refer 
to  E^0  as  a  modulating  function  in  the  time  domain.  In  the  engineering  literature,  multi¬ 
plication  by  E{0  is  called  “Quadrature  Amplitude  Shift  Keying”  or  “Quadrature  Amplitude 
Modulation”  (see,  for  example,  [26]).  “Quadrature”  refers  to  the  fact  that  the  real  and 
imaginary  parts  of  the  modulating  function  are  f  out  of  phase  with  each  other. 

Lemma  3.1  If  f  6  L2  and  /i  is  realizable ,  then 


f  *  ECoh  =  ECo  [E-bf  *  h]  ■ 

Proof  :  We  have  that 

/OO 

/(r)e2,nCo(t~T)/i(t  -  r)  dr 

-OO 

/OO 

e-2^Co Tf(r)H{t  -  r)  dr 

-OO 

=  E(.n  (£-<«/  *  /i]  ■  □ 


(20) 


Remarks  : 

1.  Thus,  to  modulate  a  channel,  we  modulate  the  incoming  signal  before  it  goes  in  and 
after  it  comes  out  of  the  channel. 

2.  The  lemma  holds  for  more  general  classes  of  signals  and  convolvers. 

3.1  Developing  a  Specific  System 

We  now  apply  Lemma  3.1  to  create  a  multichannel  system  for  a  system  with  impulse  response 

Hi{t)  = 

The  Fourier-Laplace  transform  of  /Ui(t)  is 

sin(^C) 

This  has  zeros 

*  =  {±f:teN}. 
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Let 


Then 


and 


jjL2{t)  =  Ei_n i  =  e*ltni(t) . 


MO  = 


sin(27r(C  -  |)) 

7r(C-*) 


Theorem  3.1  The  functions  pi (t) ,  n2(t)  form  a  strongly  coprime  pair  of  convolvers. 


Proof  :  We  want  to  show  that  there  exist  positive  constants  A  and  B  and  a  positive  integer 
N  such  that  i 

(y.  lft(C)|2)  ’  >  Ae-2'fl'a<l(l  +  Id)-"  ,  for  all  <  €  C. 

We  need  the  following. 

Proposition  3.1  If  19(C)!  ^  3log(2),  then  |sin(C)|  >  |e-9^. 

Proof  :  We  have  that 

|sin(C)|2  =  I  [e*>«>  +  -  2  cos  (2»(C))]  . 

If  3(C)  >  3  log(2)  >  0,  then  e29^  >  e2>  *°«(2)  =  2,  and  so 

e29(C)  +  g— 29(c)  _  2  cos(25R(C))  >2(1-  cos(2K(C)))  +  e"29^  >  e"29^  =  e-2ls^>' . 

If  S(C)  <  log(2)  <  0,  then  e-29^  >  e2^los^  =  2,  and  so 

g29(c)  g— 29(c)  _  2cos(25R(C))  >  2  (1  —  cos(25R(C)))  +  e29^  >  e29^  =  e-2!9^^  . 

Thus,  |sin(C)|2  >  |e-29^L  Taking  square  roots  proves  the  proposition.  □ 


Proposition  3.2  For  any  e  such  that  0  <  e  <  |3(C)|;  there  exists  ct>  0  such  that  |sin(C)|  > 
cee~ . 

Proof  :  We  have  that 

|sin(C)|2  =  ^  [e20(°  +  e_29(c)  -  2cos(23?(C))]  =  ^  [cosh(2$(C))  -  cos(2K(C))]  • 


2  Not  to  be  confused  with  the  zero  set  in  Section  2. 
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Now,  cosh(f)  is  continuous  and  cosh(t)  >  1  for  all  ^  €  R  with  equality  if  and  only  if  t  —  0. 
Therefore,  we  have  that  |sin(()|  is  bounded  below  by  some  constant  Ae  >  0  for  all  e  such 
that  0  <  e  <  |3(C)|-  Combining  this  observation  with  the  proposition  above,  we  get 


sin(C)|  > 


3(0  ‘f  l$(C)l  >  |log(2) 

Aee~9^  if  0  <  e  <  |9(£)|  <  \  log(2) .  □ 


Therefore,  if  |9(£)|  >  ^log(2),  |sin(27rC)|  >  fe  2lT®(0  and  sin(27r(C  —  j))  >  \e  2?r9^. 
Thus,  for  |9f(C)|  >  5log(2), 

Vin/AI2  -  1  |sin(27rC)  j2  .  |sin(27r«  -  -}))|' 

h  °  *2  ICI2  1C  -  il2 

e-47rS(C) 

>  _L  (2  lc>g(2))  +  16  -2(2ttS(C)) 

-  47r2[  (1  +  |C|)2  J 

So,  let  A  =  A-  yj (5  log(2))2  +  ^  ,5  =  1,  and  N  =  1.  Also,  for  0  <  e  <  |9(0I  <  \  log(2), 

VlO.^'H2  >  _ 4 _ 1-2(2*55(0) 

feMC)l  -  -2l(i+ici)2J 

So,  let  A  =  B  —  1,  and  N  =  1.  Finally,  for  |9(C)|  <  e,  |sin(27rC)|  >  2  |1  —  cos(4-7r5ft(C))|- 
Thus, 

5Zlft(C)|2  >  4  77  ^ >  1  \ 2  [l1  ~  cos(47rjR(C))l2  +  |1  -  cos(47t(»(C)  -  \))\2  • 
j=i  ^  LvT  +  lsIl  L  4  J. 

Now,  in  this  last  term,  the  cosines  are  —periodic  and  are  exactly  |  out  of  phase  with  each 
other.  Therefore,  this  term  is  >  (-^7=)  =  Thus, 

§ lft(c)|2  -  i  [f+W.  ' 

So,  let  A  —  B  be  arbitrary,  and  N  —  1.  Note,  all  estimates  are  independent  of  to  =  9?C- 
Combining,  we  let  B  —  N  =  1  and 

A=\t  if  0<e<  |9f(C)| 

l  f  if0<|9(0l<e. 


>  1  [1C12  +  K-!12 

-  4vr2  [  |C|2|C  -  1I2 


This  proves  the  theorem.  □ 


16 


The  next  step  is  for  us  to  construct  deconvolvers  such  that  for  an  approximate  identity 
function  ip,  we  have 

fJ>l  *  Vxrf  +  =  Ip  ■ 

We  will  need  some  computations.  We  know  that 


These  have  zeros 

Now, 

Thus,  for  z  £  Zl, 

Also,  for  z  (E  2%, 


K(C)  =  5i^>WC)  =  sE(M-J) 


ttC 


1  .  A;' 


*  =  ±2  ■*  =  4±2  ’^N' 


d  _  2cos(27rC) 

dcM c)  =  — C  'c^0- 


d 

d( 


2 

\z\ 


l/^i  (*) I  = 


1 


ttH 


Symmetric  relationships  hold  for  pi  and  its  derivative. 

Let  k  €  N,  and  let  a  =  &  +  /?,  where  0  <  /?  <  1.  We  define  the  Holder  space  Ca  of 
exponent  a  as 

{/  €  Ck :  There  exists  M  >  0  such  that  f^k\t)  —  /^(r)|  <  M  \t  —  r|^| . 

We  also  define  the  related  class  Da  as 


{/:/(«)  (1  +  M)“  €  L“}  , 
for  u  in  frequency  space.  We  have  that 


Da+l  C  Ca. 


If  /  has  compact  support,  then 

f  Gca  =>  f  e  Da. 

Choose  ip  to  be  in  the  Hdlder  space  C3+v,  r)  >  0,  with  support  in  (—2,  2)  such  that  ip  >  0  and 
ip(t)  dt  =  1.  We  refer  to  ip  as  the  auxiliary  function  of  the  construction.  The  conditions 
that  we  impose  on  ip  make  it  an  approximate  identity  function.  Its  smoothness  and  the  size 
of  its  support  guarantee  that  the  deconvolvers  are  compactly  supported.  For  the  compactly 
supported  function  ip,  we  have  that  if  rj  >  0  and  ip  €  C3+n,  then 

l*wl  5  (i  + 
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for  z  G  Zi  U  Z2. 

We  will  assume,  for  technical  reasons,  that  /  G  C°°  n  L2(K).  We  will  use  a  density 
argument  at  the  end  of  this  section  to  extend  the  result  to  all  of  L2.  Let 

(fA)  =  [  f(t)i/>(t)dt 

JR 

be  the  L 2  inner  product  of  /  and  ip,  and  let 

{v;  i>)  =  /  i>{t)  v{t)  dt 

JR 

denote  the  dual  product  between  ip  G  £  and  v  G 

We  want  to  show  that  there  exist  a  sequence  of  concentric  circles  Tn  with  centers  at  the 
origin  and  radii  pn  such  that  limn^.00  pn  =  oo  and  such  that  there  exists  c  >  0  with 

|sin(27rC)|  >  ce2lT^ 


and 

sin(2-7r(C  -  ^))  >  ce2*^ 

for  all  C  G  rn.  The  choice  of  such  Tn  is  obvious.  Let 

f  .  .  2n  +  1 
=  |c :  |C|  =  — g —  >  n  £  N 


(This  is  exactly  half-way  between  the  zeros  of  p\  and  p2.)  For  (  G  Tn,  we  have  that 


f  if  l»«)l  >  Uog(2) 

2  if  i^(c)i  <  I  iog(2> 


We  choose  c  =  ^,  which  is  independent  of  n. 
Theorem  3.2  The  set 


Mi (t)  =  X[-i,i](<)  ,  p2(t)  =  EiXhhl](t) 
is  a  strongly  coprime  pair  of  convolvers.  The  deconvolvers  such  that 


(f,  V)  =  {"hri  f  *  Mi)  +  /  *  M2) 


are  given  by  the  formulae 


VlAt)  - 
*At) = S, 


(21) 


(22) 
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Also, 


0(0  — >  <5  as  supp(0)  — >■  {0} 

and 

•0(C)  — >  1  as  supp(0)  — >  {0} 

in  the  sense  of  distributions. 

Remark  :  The  formulae  for  the  are  a  type  of  Lagrange  interpolation,  where  the  analytic 
functions  0(C)  were  constructed  from  known  values  on  a  discrete  set  of  data.  The  Cauchy 
residue  theory  is  used  to  convert  these  discrete  values  to  analytic  functions,  in  the  same 
fashion  as  sampling  formulae  convert  discrete  functions  to  analog  functions.  The  formulae 
are  also  related  to  classical  Shannon  sampling,  with  the  functions 

i  (MOV 
\C  -  z) 


acting  as  interpolators. 
Proof  :  We  have  that 


MO  = 


sin(27rC) 

ttC 


,  MO  = 


sin(27r(C  -  Q) 

AC  -  |) 


dC 


Pi{z)  ^0, 


and 


Mz) 


2 

\z\ 


for  zeZi.  ^ 

Since  0(C)  is  entire,  we  may  represent  0  by  the  Cauchy  integral  formula,  i.e., 


Lemma  3.2  For  C  €  {C  :  |C|  <  Pn), 

0(C)  =  MO  E 


zSZ2 

\*\<Pn 


+  MO  E 


z€Z  ! 
1*1  <Pn 


i  9(z)dz. 
Jrn  z  -  C 

0(0  j 

(  MO 

MO^MO 

{C-z 

0(0  j 

(MO 

MO^MO 

U-* 

+  MQ, 

where  Rn(Q  — »  0  as  n  — »  oo. 
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Proof : 


$(C) 


J-i  it 

2ni  Jvn  z  — 


dz 


—  / 
2tt  i  Jr 


$(*)_ 

C 

$(z)fii(z)fflz)  -  $(z)JIi(C)MC) 


+ 


(z  -  0  Mz)H2 {z)) 


±f 

2iti  Jrn 


rn  {z  -  C)  (ni{z)n2(z)) 

${z)  \pi{z)fe(z)  -  MOMQ]  dz 


dz 

dz 


(z-Q  (Mi(*)M2 {z)) 

+  ^M)MO  £  (z  _  {)  (^(i)|r2(z)) iz- 

Now,  for  z  €  Zi  or  z  €  Z2,  (Hi(z)  ^2(z))  =  0,  but  ^  (hi(z)/j,2(z))  7^  0.  Thus,  the  function 


( 


)  has  simple  poles  at  z  G  Z\  U  Z2.  Therefore,  has  simple 


Hi(z)m(z) 

poles  in  Zi  U  Z2 ,  and  so  by  the  Cauchy  Residue  Theorem, 


/ 

•/r„ 


$(*)  -  MCWQ] 


z€22 

|*1<P« 


!*!<<>■ 


(«-C)  (Mi {z)to{z)) 

VKf) 

A 
dC 

i>{z) 

A 

dC 


dz 


^l(C)  ^  Hx{z)jzHi{z)  [C-zJ 


(fS) 

(M0\ 


+  ^  ,5,  \C~z)' 


Let 


*«>  -  5Sfi«>«<«  j£.(TT 


$(*) 


dz. 


(*-  0  (M1WM2 (*)) 

To  finish,  we  need  to  show  that  /^(C)  — >  0  as  n  — >  00.  We  were  able  to  choose  Tn  such 
that  for  z  €  Tn, 

|sin(27rz)|  >  ^e2^', 


sin(27r(z  -  -)) 


>/2 


,27r|3z| 


By  the  Paley-Wiener-Schwartz  Theorem,  there  exist  positive  constants  C  and  A ,  with  A  <2 
and  C  =  C(N),  such  that  for  any  positive  integer  N  and  any  z  G  C, 


${z) |  <C{l  +  \z\)-N<?"AW 
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Therefore,  for  £  fixed  and  any  fixed  iV  >  0, 


\Rn(0 1  < 

< 

< 


MOMOI 


1 

I 

■0(0 

2m 

Jrn 

(z  —  0  (M0M0) 

1 


MOMOI  sup 

zer„  2n 


M) 


MOMOI  sup 


(z  —  0  (/A  (0/^2 (0) 

C(1  +  \z\)~N  pne^A^z\ 


dz 

2npn 


MOMOI 


(IM-ICH)  (i^l«»l) 

+  Pn)~N Pn  2t x{A- 2)pn 


|(Pn-|C|) 


Since  A  <  2,  this  last  quantity  — >  0  as  n  — >  00.  This  completes  the  proof  of  the  lemma. 

□ 


To  finish  the  proof  of  the  theorem,  we  need  to  show  that  the  series  expressions  in  the 
previous  lemma  converge  in  the  sense  of  £',  that  is,  as  analytic  functions  which  satisfy  the 
PITS'  growth  estimate. 


Lemma  3.3  The  series 


y'  0(0  f  MOO 
kii  MO^MO  ~  z' 

1*1  <Pn 


i,j  =  1,2 ,  i^j 


converge  in  £' . 


Proof  :  If  MO  =  0,  then 


Also,  if  MO  =  0, 


*  =  1,2. 


MO  I 


1 

7T  \z 


Mi- 


The  functions  MO  are  entire.  Therefore,  on  every  compact  subset  K  of  C,  there  exists 
a  constant  c(K)  such  that  |M0I  <  c(M  f°r  C  £  K.  The  function  ip( t )  is  C3+n  for  17  >  0, 
and  therefore  its  transform  will  decay  on  the  real  axis.  In  particular,  for  z  £  Zi  L)  Z2,  by 
PITS  there  exists  a  constant  d  —  d( 3  +  rj)  such  that 


0(0 1  < 


d 


(1  + 10) 


3+t)  • 
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Therefore, 


< 

< 


R*)|  1/5(01 

\MZ)^MZ)\  |C“*| 

nd  2  c(K) 

2 (i  +  \z\)3+T}  z  MCI- Nl 

c(K)-Kd  1  .  .  . 

2\z\"  ||d -Ml  *’J  ’  ’  '±1 


In  particular,  for  (  fixed, 


$(z) 

1/5(01 

MZ)«MZ) 

\C-z\ 

< 


(constant)^ 


fi+v 


Thus,  by  the  Weierstrass  M-test,  the  series 


fe)  \m)\ , 


S 


*,i  =  1,2  ,  i  ^  j 


\z\<Pn 


both  converge  uniformly  on  compact  subsets  of  C,  and  therefore  represent  entire  functions. 
Now,  for  each  z  £  and  |(  -  z\  >  e  >  0  with  e  <  1/2,  there  exists  a  constant  a{ 


MO 


If  |C  —  z\  <  e,  then 


C-z 

MO 


<ilffi(OI<-«2’"acl- 


\C-z 

Therefore,  there  exists  k,  >  0  such  that 

MO 


<  sup 


£m) 


C-z 


<  Kilt*™ 


independently  of  z  €  Z{. 

From  the  estimates  above,  we  have  that  the  series 


£  \$(z)\ 

**i  \pj(z)ffi(z)\ 
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is  bounded  by  some  B.  Therefore 


^(z) _ Mt(C) 

{ C~z) 

\A<Pn 

<  y'  IV’Wl  1/5(01 

<  BKie2^c  I, 

Thus,  the  series  converges  in  £' .  This  completes  the  proof  of  the  lemma  and  the  theorem. 

□ 

We  now  give  explicit  formulae  for  the  deconvolvers  without  using  inverse  transforms. 


Lemma  3.4  If  p,((3)  =  0,  then 


«-« 


=  no. 


where 


(T(/?);^)  =  "2 *i  J  (J*  n{x)dx. 


Proof  :  Since  T  €  £', 


f  =  (T]e~27rit<:) 

—  —2ni  £(f‘-  e27rz/3(t  da; 

=  -2tt*  j°°  ix(x)e-2™Pdx 

/oo  /  p-27rixC„2OTt®/?  _  1  \ 


“  L  I"2"  — 2irt  (C  —  0) 

/•°°  n(x)e~2irixi  j  f°°  n(x)e~2mxP  ^ 
=  J-  oo  (C-0)  J-oo  (C-P) 


J-co  (C  —  /?) 
js(0  -  Hg) 
«-/?) 
m 
(c-py 


since  /i(/?)  =  0.  □ 


Lemma  3.5  For  Hi  as  above  and  z  €  i,j  =  1,  2, 


fc(C) 

(C-*) 


=  35(0, 
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where 


m)  = 


T,(i)  =  -  (e2’“<‘+‘)  -  l)  %(_!,,](*) , 

z 

2*i(z-\)(t+l)  _  A  _ 


(-1) 


Proof  :  For  z  €.  Zx, 

{Tx\ip)  =  -27ri  J  Qf  ip(t)e2irlz{t~x)  dt)  X[_ltXfx)dx 
=  -2m  J1  (jT  ip(t)e2*iz{t~x)  dt)  dx 

=  2m  J°  (ip(t)  J*  e2™V~x)  dx)  dt 
+ 


-  2m  jf 1  (i>{t)  jf  *  e2niz(-t~x)  dx)  dt 
I  (e2«,(t+i)  _  !)  X[_M ](t);^\ 


For  z  €  i?2, 


<T2;V>)  =  -27ri|°°  ^^(t)e2,rfo(t-a!)(ft)e(7r/2)feX[_1,1](*)da: 

=  -2ttz  y)  (jT  ip(t)e2”iz{t-x)  dt)  e^/2)ix  dx 
=  2m  f  (ip(t)e2nizt  y)  e"27ri((*"4>a;)  dx)  dt 

+  -  2m  y1  (V^e2™*  f  e-27ri((*-*)x)  dx)  dt 
=  ((JTjJ  (e27riM)(m)  -  l)  □ 


These  last  two  lemmas  give  us  explicit  formulae  for  the  i/j  in  the  time  domain.  We  add 
that  the  estimates  given  in  the  proof  of  the  previous  theorem  give  us  that  ip  only  needs  to 
be  a  compactly  supported  C3+v  function  for  r)  >  0,  because  this  is  sufficient  to  guarantee 
that  the  series  representations  for  the  ut  converge  to  compactly  supported  functions. 

Theorem  3.3  Let  f  €  C'00nL2(R);  and  for  r)  >  0  let  ip  be  a  function  with  support  in  (-2, 2) 
in  the  Holder  space  C'3+?J  such  that  ip  >  0  and  ff0^  ip(t)  dt  =  1.  Given 

Ah (t)  =  X[_i,i](t)  ,  A toft)  =  e*ltX[_iti ](t), 
the  deconvolvers  Ui ^  such  that 

f*ip  =  (f*Hi)*  vx^  +  (/  *  Ah)  *  v^> 
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are  given  by  the  formulae 


*At)  =  E  W=-( , 
•%*(*)  =  E  , 


(23) 

(24) 


77ie  function  f  *  ip  is  an  arbitrarily  close  approximation  of  f  which  converges  to  f  in  the 
sense  of  distributions  as  supp(^)  — »  {0}. 


To  finish  this  section,  we  need  to  show  that  the  deconvolution  procedure  given  in  Theorem 
3.3  works  for  general  L 2  functions.  We  first  note  that  C'0OnL2(M)  is  dense  in  L 2 .  We  also  note 
that  the  functions  are  compactly  supported  and  continuous.  Therefore,  since  f*pt  G  T2, 
f  *  Pi*  ^  G  L 2  with 

\\f*Pi*^h  <  ||/*A*i||2||*/<^||i- 

Similarly, 

ll/*^lla  <  ll/IMMli- 

Define  the  operator 

V:  L2(R)  — >  L2(R) 
by 

T>(f)  =  f*1p~{f*pl*  Virf  +  f*P  2  *  ^2, ip)  • 

The  operator  V  is  clearly  linear.  The  estimates  above  give  us  that  V  is  bounded,  and 
therefore  continuous.  Since  V  —  0  on  C°°nL2(R)  and  is  continuous,  V  —  0  on  L2(R).  Thus, 
the  deconvolution  procedure  works  for  general  functions  of  finite  energy. 

Remark  :  Thus,  given  a  fixed  system,  it  is  possible  to  modify  the  system  via  some  easily 
performed  transform  to  create  a  strongly  coprime  system,  and  consequently  recover  the 
complete  input  function.  This  modulation  technique  works  for  creating  strongly  coprime 
systems  for  B-spline  systems.  In  a  similar  vein,  given  a  strongly  coprime  system 

is  created  by  rotating  the  square  by  j  and  j  (see  [13]). 

Remark  :  The  deconvolvers  are  not  unique.  For  example,  given  two  deconvolvers  and 
u2,  A  G  R,  and  a  compactly  supported  function  p,  the  pair 

+  \r)  *  p2  ,  i>2  —  -  \r)  *  pi 

is  also  a  set  of  deconvolvers.  The  non-uniqueness  of  the  deconvolvers  allows  one  to  de¬ 
velop  deconvolvers  that  are  optimal  with  respect  to  a  given  condition.  In  [13],  deconvolvers 
that  optimize  the  signal  to  noise  ratio  of  the  signal  relative  to  white  Gaussian  noise  were 
constructed.  Other  types  of  noise  could  be  dealt  with  in  similar  fashion. 

Remark  :  We  can  also  use  a  family  of  zero-mean  Gaussians  for  our  auxiliary  function  ip.  If 

i’xif)  =  ^=e_t2/(2A2)  • 
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Figure  3.  Gaussian  Function  G  Approximate  Deconvolver  System. 

3. a.)  Channel  1  Deconvolver  v\. 

3.b.)  Channel  2  Deconvolver  v2. 

3.c.)  G  *  *  V\. 

3.d.)  G*  n 2*  v2. 

3.e.)  G  —  G  *  ^  *  V\  +  G  *  ^2  *  ^2,  the  resulting  approximate  delta. 
3.f.)  The  Fourier  transform,  or  frequency  response,  of  this  system. 


Figure  4.  Transfer  Response  for  the  System. 

4. a.)  The  response  of  each  channel,  in  the  two  channel  system  of  Figure  3,  to  two  adjacent 
Gaussian  pulses  (denoted  as  input  signal  /). 

(fi  1  *  f  —  solid  ,  H2*  f  ~  dotted) 


4.b.)  Transfer  response  for  the  system  given  in  Figure  3.  Here  the  dotted  line  represents  the 
input  function  while  the  solid  line  shows  the  system  response. 
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4a.  The  response  of  each  channel,  in  the  two  channel  system  of  Fig.  3,  to  two  Gaussian  pulses 
(denoted  as  the  input  signal,  f). 
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4b.  The  system  transfer  response,  for  the  two  channel  system  of  Fig.  1 ,  when  the  input  function  is 
the  sum  of  two  time  shifted  Gaussian  pulses.  The  dotted  line  in  the  input  function,  the  solid  line  in 
the  deconvolved  function. 


Figure  4.  Transfer  response  for  the  system  of  Figure  3. 


then 


iMC)  =  ■ 

We  restrict  A  >  4,  and  so  twice  the  standard  deviation  is  <  2.  As  A  — >  oo  ip\  — »  1  in  the 
sense  of  convergence  in  S'  (ip\  is  the  transform  of  an  approximate  identity) . 

Convergence  estimates  must  now  be  in  terms  of  S'  and  S'.  The  deconvolvers  developed 
by  using  are  not  compactly  supported.  However,  from  a  numerical  viewpoint,  they 
are  “essentially  compactly  supported,”  in  that  for  t  >  4A,  the  numerical  values  of  the  de¬ 
convolvers  essentially  equal  zero.  Figures  3  and  4  repeat  the  simulation  of  Figures  1  and  2 
using  a  Gaussian  auxiliary  function. 


3.2  Sampling  Revisited 

We  first  need  to  recall  the  following  theorem. 

Theorem  3.4  ([32])  Let  F  be  a  function  of  finite  energy  on  [— T,  T]  (F  €  L2[—T,  T ]).  Then 

m  =  A  E  F(n/2T)e2"<’*'27')'X|-T,T|W . 

Moreover,  given  any  sequence  {an}nez  C  C  with  |on|2  <  oo  (i.e.,  {an}  6  l2),  then  there 
exists  a  function  F  €  L2[—T,T ]  such  that 

F(n/2T)  =  an 


for  all  n  £  Z.  In  fact, 


Fit)  =  E  F(n/2T)e^/^‘XhT,rl  (t) 


neZ 


Also, 


£  F(t  +  2nT)  =  ft  £  F(n/2T)e2'i<“/2T)‘X|_T,T|(() 
n€Z  neZ 

where  the  sums  converge  in  L2[—T,T\. 


Let  k  €  N.  The  zeros  of  /ii(C)  =  sm[2^  are  Z\  =  {if}’  wh*le  f°r  A*2 (t)  =  Eifj, i 


e?rtyUi(*)>  Afe(C)  =  Sml2ffl)4))  and  so  Z2  =  {|  ±  |}.  Thus,  for 


n  G  Z, 


Now,  the  construction  of  the  deconvolvers  in  the  previous  section  assumed  that  the  auxiliary 
function  % b  had  support  C  (-2,2).  A  construction  of  ip  via  Shannon  sampling  in  frequency 
is  possible.  The  Nyquist  rate  is  \  -  exactly  the  same  rate  as  Z\  (J  Z2.  These  observations 
yield  the  following. 
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Theorem  3.5  Let  f  £  L2(R),  and  for  r)  >  0  let  ip  be  a  function  with  support  in  (—2, 2)  in 
the  Holder  space  C3+v  such  that  ip  >  0  and  ip(t)  dt  =  1.  Given 

Mi (t)  =  X[_x>1](t)  ,  p2(t)  =  e^X^qit), 
the  deconvolvers  such  that 

/  *  V*  =  (/  *  Mi)  *  vi,*  +  (/  *  M2)  * 


are  given  by  the  formulae 


"lA*)  =  \ 


kx  +  J2  W2)l  girint 


n^O 


Ml  ((l/4)  +  (n/2)) 


]{t) , 


(25) 


"2M  =  g 


M2(n/2) 


^[-1,1]  CO ) 


(26) 


where 


#((1/4)  +  (n/2)) 
/4)  +  (r 

$(«/2) 


S(  ir+W/4)  +  (n/2))  ’ 


*-£(-irW 

The  function  f  *  ip  is  an  arbitrarily  close  approximation  of  f  which  converges  to  f  in  the 
sense  of  distributions  as  supp (ip)  — >  {0}. 

Proof  :  By  Theorem  3.1,  and  p2  form  a  strongly  coprime  pair.  Therefore,  there  exists  a 
positive  constant  C  and  a  positive  integer  N  such  that  for  all  u>  €  R, 

1  Mini  +  imiMi  >  c (1  +  m)  N . 

Thus,  since  fi{  =  0  and  M2  (|  +  2)  =  0  for  all  n  €  Z  \  {0}, 


Mi  +  -)  >C(l  +  |n|)  N  , 


n 

M2(2 


>  C  (1  +  |n|) 


-N 


for  all  n  e  Z  \  {0}.  Moreover,  estimates  developed  in  the  proof  of  Lemma  3.3  give  us  that 
N  =  2.  Now,  since  ip  is  a  function  in  the  Holder  space  C3+v, 


\$(z)\  < 


d 


(i  +  N) 


3+J) 
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for  z  G  Z\  U  >^2,  and  so 

ft(l/4)  +  (n/2))  v  #((1/4)  +  (n/2))  2 

^((V4)  +  (n/2))  <  ’4-  ^((1/4)  + (n/2))  <0°’ 

#(n/2)  y'  #(n/2)  , 

Hence,  the  sums  in  the  definitions  of  z/^,  z^  converge  absolutely  and  uniformly,  and  the 
coefficients  of  the  deconvolvers  are  in  £2. 

Now,  for  z  £  Zx  U  22, 

fii (z)vQ>(z)  +  Mz)^(z)  =  $(z)  >  (27) 

for  if  z  £  Z{ ,  /^(z)  =  0  and  Vj^j(z)  =  4^,  for  i  j,  i,j  =  1, 2.  Also,  for  z  =  0  and  z  =  1/4, 
the  formulae  agree.  Thus 

tH(n/4)vQ,(n/4)  +  /Za(n/4)i^(n/4)  =  #(n/4)  (28) 


for  all  n  €  Z.  Letting 

F  =  + 

we  have  that  F  is  continuous  and  satisfies  the  hypotheses  of  Theorem  3.4.  Therefore,  F  =  0, 
which  proves  the  result.  □ 

We  can  again  reconcile  the  apparent  differences  between  the  deconvolver  formulae  by 
computation.  We  have  for  z  €  Zi 


.  2 cos  (2tt  (z  —  |)) 

AC  ^2(2)  —  /  i\ 

d(  [z  -  J) 


Thus, 


"lAt) 

v-  #(*)  1 


=  V  Ti — (  7  1  \  (e27ri(*~<)(H1)  _  i')  e("/2)«x  ,(i)  ] 

=  V  (  7  1  v  (V™M)«+D  -  l)  e^/^Xr. 

Ah(*)  2  cos  (2tt  (z  -  J))  \  (*  -  J)  V  '  ’  / 


zez2  l*nz)  2  cos  (27T  (z  - 
1  „  m/4)  +  (n/2)) 


-  2  §a«i/4)  +  (n/2))(~ir  (e2’i<"/2><‘+1)  -  ')  eW2)“"'-'’ll(t) 

=  3  E  Smm 

=  5  £  ((-1}"+1 + e’”‘)  • 
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Similarly,  for  z  G  Z\ 


Thus, 


2  cos(27tz) 
z 


v*  At) 

V'  1  /'„27rw(t+l) 

M*)  fcM*)  ^ 

_  $(Z)  £  /"l  /e2ir**(i+l) 

z^Mz)2cos(z)\z^ 


1)  X[-i,n«) 
l)  X,-!, ,](<)) 


1  ^  ^(n/2) 

2  p^(n/2) 

1  ^  $(n/2) 

2  /£(n/2) 

1  ^  #(n/2) 

2  ^(n/2) 


( — l)n  (e27ri(n/2)(t+1)  -  l)  X[_M](t) 
(  l)n  ((e7r*)ne7rint  -  l) 

((-l)n+1  +  e™‘)  XH;1](t). 


Remark  :  The  development  of  the  deconvolvers  via  the  Jacobi  interpolation  formula  and 
Cauchy  residue  theory  is  now  a  well-developed  tool,  given  the  theoretical  base  established 
by  Berenstein,  Gay,  Taylor,  Yger,  et  al.  It  has  a  flexibility  in  the  one-variable  case  that 
allows  for  its  use  in  not  only  general  deconvolution  problems,  but  also  in  the  development  of 
filters,  etc.  The  key  to  this  is  the  flexibility  of  the  Cauchy  residue  calculus  in  one- variable. 
As  is  well-known,  the  story  in  several  variables  is  different.  Berenstein,  Gay,  Taylor,  Yger, 
et  al.  have  given  us  working  formulae  for  specific  situations.  There  are  no  general  formulae 
in  several  variables  for  the  computation  of  the  needed  residues. 

Sampling  lies  juxtaposition  these  methods.  It  does  not  have  the  tremendous  flexibility 
of  the  complex  methods  in  one  variable,  but  it  also  does  not  carry  with  it  comparable 
computational  difficulties  in  several  variables. 

Both  methods  will  be  used  in  further  developments  of  the  theory. 


4  Discussion  of  Applications 

We  now  give  a  brief  discussion  of  the  applications  of  these  deconvolution  techniques,  and 
refer  to  [13],  [14],  and  [15].  We  also  mention  [17].  The  multichannel  theory  was  extended 
to  a  much  broader  class  of  convolvers  in  that  paper,  and  was  used  to  design  filters  in  some 
specific  cases. 

Deconvolution  is  applicable  in  any  area  of  signal  and  image  processing  in  which  it  is 
desirable  to  know  a  high  degree  of  detail  from  the  input  data.  A  multichannel  system 
consisting  of  an  array  of  strongly  coprime  convolvers  {ni}  guarantees  that  any  information 
about  the  input  function  /  lost  by  one  convolver  is  retained  by  another  in  the  array  as 
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Si  —  f  *  Hi  for  some  i.  The  signals  Si  are  then  filtered  by  the  and  added,  resulting  in 
the  reconstruction  of  /.  Again,  the  methods  are  both  linear  (convolution  with  deconvolvers) 
and  realizable  (the  support  of  the  deconvolvers  being  contained  in  the  bounded  support 
of  the  convolvers),  with  deconvolution  at  a  point  t  G  Rn  depending  only  on  data  near  t. 
The  theory  assumes  no  a  priori  information  about  the  input  signals.  The  non-uniqueness 
of  the  deconvolvers  allows  one  to  develop  deconvolvers  that  are  optimal  with  respect  to  a 
given  condition.  We  also  have  precise  knowledge  about  digital  versions  of  the  deconvolvers, 
although  the  sampling  rate  over  those  domains  still  depends  on  the  bandwidth  of  the  received 
signals,  which  equal  the  bandwidth  of  the  original  signal.  Additionally,  we  can  combine  this 
linear  procedure  with  other  linear  procedures  which  extract  information  about  the  signal, 
e.g.,  discrete  wavelet  transforms  (see  [16]). 

4.1  Image  Processing 

Deconvolution  will  be  very  useful  in  image  processing.  In  a  system  designed  to  perform 
image  analysis,  the  deconvolvers  would  act  as  an  initial  enhancing  filter  for  the  image.  After 
this  process,  algorithms  such  as  edge  detection,  etc.  would  then  be  working  with  data  arbi¬ 
trarily  close  to  the  original  data,  instead  of  data  in  which  the  high  frequency  information 
has  been  lost.  Further,  these  various  types  of  processing  could  possibly  be  combined  with 
deconvolution  so  as  to  produce  the  results  in  a  single  processing  step. 

The  problem  of  getting  an  exact  representation  of  a  pixel  image  is  addressed  in  [13]  and 
[14].  Here,  an  image  is  restored  from  the  data  gathered  by  a  set  of  three  photo-detectors. 
Each  photo-detector  is  modeled  as  the  integrator  over  time  of  an  image  /  in  a  compact  region 
of  R2.  Thus,  the  model  for  the  image  data  is  s  =  /*//,  where  p  is  the  characteristic  function 
X  of  some  bounded  region  in  M2.  The  strongly  coprime  system  is  pi  =  = 

X[-V2,V2]x[-V2,V2]  >  ^3  = 

To  deconvolve,  it  is  necessary  to  multiplex  in  time,  i.e.,  the  time  interval  used  for  inte¬ 
gration  will  be  subdivided  and  used  to  obtain  finer  spatial  sampling  in  multiple  channels. 
Although  no  larger  time  interval  is  required,  the  use  of  the  time  interval  for  multiplex¬ 
ing  means  that  we  are  not  deconvolving  the  time  integration  but  rather  depending  upon 
a  relative  time  invariance  of  the  signal  over  this  interval.  Two  levels  of  multiplexing  are 
required  for  detector  arrays:  one  multiplexing  is  to  approximate  a  convolution,  the  second  is 
to  obtain  multiple  convolutions.  The  deconvolvers  used  are  two  dimensional  versions  of  the 
deconvolvers  described  in  Section  2.  These  deconvolvers  are  separable. 

For  these  detector  applications  an  important  topic  both  for  hardware  development  and  for 
models  of  vision  is  the  design  of  multiplexed  arrays.  We  want  to  know  how  to  best  arrange 
in  a  single  array  a  set  of  approximately  strongly  coprime  detectors,  how  to  best  “scan,” 
“dither,”  or  “sweep”  the  array  in  order  to  obtain  the  desired  oversamples,  and  how  to  use 
these  measurements  for  deconvolution  and  signal  estimation  as  well  as  for  the  determination 
of  information  about  the  scene  (source  of  the  signal)  from  the  signal.  Associated  with  these 
questions  are  the  obvious  ones  of  knowing  or  estimating  temporal  and  spatial  phase  of  sample 
points  in  the  array.  Because  of  such  considerations,  the  strongly  coprime  pair  of  convolvers 
given  by  the  characteristic  functions  of  two  disks  with  suitably  different  radii  have  nice 
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properties  as  well  as  obvious  analogs  in  biological  retinas.  The  strongly  coprime  condition 
is  satisfied  with  only  two  disks  of  radii  ri,  7*2,  with  r\fr<i  =  k  G  N.  Moreover,  two  such  disks 
can  be  positioned  concentrically,  which  eliminates  one  phasing  problem.  In  addition,  this 
configuration  may  reflect  the  arrangements  in  biological  retinas  that  are  described  by  Marr 
[27]  as  the  “central  surround”  response  used  (presumably)  for  edge  detection. 

The  results  on  deconvolution  are  for  distributions  with  compact  support.  The  transforms 
of  these  distributions  are  entire  functions  and  hence  have  at  most  isolated  zeros.  Further, 
these  results  are  about  a  set  of  these  distributions  such  that  no  isolated  zero  is  common  to 
all  in  the  set.  Obviously,  these  deconvolution  results  are  not  applicable  to  convolvers  b  that 
are  band-limited.  They  are,  however,  applicable  to  the  compactly  supported  convolvers  in 
a  system:  if  the  yU;,  i  =  1, . . .  ,  n,  are  strongly  coprime  and  each  /i*  is  the  composition  of  all 
of  the  compactly  supported  convolvers  of  the  ith  linear  system,  and  if  b  is  the  composition 
of  the  band-limited  convolvers  for  each  system,  then,  with  deconvolvers  i  =  1,...  ,  n, 
Ya=\  b  *  fa  *  Vi  =  b.  That  is,  we  can  reconstruct  up  to  the  band-limited  response. 

A  primary  class  of  sensor  components  that  are  modeled  by  band-limited  convolvers  is  the 
class  described  by  an  aperture  response  function  for  electromagnetic  or  acoustic  waves,  such 
as  the  point  spread  function  of  an  optical  lens.  (The  transform  rather  than  the  convolver 
has  compact  support.)  The  band-limit  is  the  diffraction  limit.  To  apply  the  theory,  we 
would  want  compactly  supported  strongly  coprime  convolvers  in  a  system,  and  construct  up 
to  the  band-limit.  An  example  is  given  by  an  optical  imaging  system  for  incoherent  light 
(e.g.,  a  telescope).  The  system  can  be  modeled  as  the  convolution  7*  =  h  *  Ig  [23]  ,  where 
I  is  intensity,  and  the  impulse  response  is  determined  by  the  pupil  function  of  the  system 
(the  modulus  squared  of  the  Fourier  transform  of  the  pupil  function).  In  many  cases,  the 
impulse  response  is  radially  symmetric  and  given  by  Bessel  functions.  For  these  systems,  the 
techniques  used  in  solving  the  Pompeiu  problem  will  have  application. 

To  extend  the  knowledge  of  the  image  past  the  band-limit  requires  both  an  assumption 
about  the  image  and  some  other  techniques  in  reconstruction.  For  systems  in  which  the 
impulse  response  is  not  compactly  supported,  to  get  information  from  the  system  at  any 
given  point,  one  would  have  to  assume  that  most  of  the  information  contained  in  the  image 
is  contained  in  a  compact  region.  (Otherwise,  to  know  the  system  output  anywhere  requires 
integration  over  all  of  R2.)  If  we  assume  that  the  image  is  compactly  supported  (certainly 
a  reasonable  assumption  for  most  imaging  systems),  the  transform  of  the  image  or  image 
intensity  is  an  entire  function  in  two  complex  variables.  The  transform  of  the  system  gives 
us  that  we  know  this  analytic  function  inside  some  ball,  determined  by  the  system  band¬ 
width.  Moreover,  the  function  is  corrupted  by  noise.  To  extend  our  knowledge  of  the  image 
past  the  system  bandlimit,  we  could  in  theory  analytically  continue  to  a  larger  disk,  and 
transform  back.  We  are  examining  the  following  technique  for  reconstructing  the  function, 
which  involves  using  the  Bergman  reproducing  kernel.  Since  the  image  is  L 2  and  compactly 
supported,  its  transform  is  L2  and  entire.  The  transform  can  therefore  be  expanded  in  terms 
of  the  kernel  -  if  t  =  u  +  iv,  f(z)  =  f  f  f(t)K(z,t)dudv.  The  Bergman  basis  functions  cun  for 
a  ball  are  known  explicitly.  We  can  expand  the  original  function  in  a  ball  larger  than  that 
determined  by  the  bandwidth.  Computation  of  the  Bergman  representation  is  approximated 
by  expanding  /  in  terms  of  these  basis  functions,  i.e.,  /  «  In  this  com- 
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putation,  the  integration  will  have  the  effect  of  averaging  the  noise.  Exact  estimates  need  to 
be  computed. 


4.2  Signal  Processing 


Multichannel  deconvolution  theory  will  be  useful  in  radar  and  sonar.  The  sensor  components 
which  are  modeled  by  compactly  supported  convolvers  include  nearly  all  of  the  various  pulses 
used  in  radar  and  sonar.  The  convolvers  describe  the  modulation  of  the  transmitted  carrier 
frequency  and  the  deconvolvers  describe  the  filters  used  to  process  the  demodulated  return 
signal.  The  deconvolution  theorem  says  in  this  context  that  if  a  set  of  strongly  coprime 
pulses  are  transmitted  and  if  each  pulse  of  the  set  is  separately  processed,  then  the  strongly 
coprime  pulses  can  be  completely  compressed,  without  sidelobes,  in  the  case  of  zero  (or 
known)  doppler  shift.  The  types  of  pulses  that  can  be  used  to  make  up  a  strongly  coprime 
set  include  any  of  the  pulse  types  used  in  radar  and  sonar.  Pairs  of  strongly  coprime  chirps 
is  one  such  example.  This  feature  suggests  that  it  is  possible  to  consider  the  enhancement 
of  the  range  resolution  performance  of  pulse  types  that  perform  optimally  relative  to  some 
other  performance  measure. 

The  discussion  of  band-limited  convolvers  given  above  is  also  applicable.  The  beam 
pattern  for  acoustic  transducers  and  for  radar  antennas,  as  well  as  for  the  real  and  synthetic 
arrays  of  either,  is  a  convolver  whose  Fourier  transform  is  band-limited.  To  apply  the  theory, 
one  would  want  compactly  supported  strongly  coprime  convolvers  in  a  system. 

We  can  apply  the  deconvolution  theory  to  the  acoustic  signals  f(x,y,z,t).  The  output 
signals  Si(x,  y,  z,  t )  =  f(x,  y ,  z,  t)*Hi(x,  y,  z)  vary  with  both  the  position  of  the  sensors  and  the 
time  at  which  the  signal  is  received.  The  theory  has  to  be  suitably  modified  to  take  this  prob¬ 
lem  into  account.  Ideally,  we  would  like  to  have  a  Green’s  function  G(x,y,z,t]X0,yo,Zo,t0) 
for  our  model  of  the  atmospheric  effects.  Then,  f(x,y,z,t)  is  represented  by  an  integral 
equation  with  terms  G,  the  initial  signal  information  g  and  | \  =  h,  atmospheric  inhomo¬ 
geneity  A,  and  boundary  conditions  -  namely,  f  =  ^*g  +  G*h  +  JtG*  A  +  boundary  terms. 
If  the  boundary  effects  can  be  repressed,  then 


f(x,  y,  z,  t)  =  dG("X^  Z ’  —  *  g(x,  y,  z) 


dt 


+ 


G(x,y,z,t)*h(x,y,z)  +  J  G{x,y,z,r)  *  A(x,y,z,r). 


If  /  is  overdetermined  by  a  strongly  coprime  set  of  sensors  {//;}  as  {s^  =  f  *  //,},  convolution 
with  the  sensors  is  over  the  spatial  variables  a;,  y,  z.  Then  for  as  above,  and  as  convolution 
is  over  x ,  y,  z, 


i= 1 


"  (dG  „  , 

=  2_,  *  g  +  G  *  h  +  J  G  *  AJ  *  fa*  Vi 

=  *  g  +  G  *  h  +  J  G  *  A^j  *  ui) 

=  ^ *f  +  G*h+[tG*A . 
at  J 
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Finally,  in  [17],  we  explored  the  idea  that  by  using  duality  and  by  solving  a  modified 
Bezout  equation,  we  may  use  filters  in  one  class  to  create  a  filter  in  another,  for  example, 
to  create  an  ideal  low-pass  filter  out  of  truncated  sine  filters.  If  //*  =  sincaj,  for 

ai  =  1,  a2  =  y/p,  p  prime,  then  {/q,//2}  form  a  strongly  coprime  pair.  Therefore,  in  the 
time  domain,  we  can  solve  the  modified  Bezout  equation 


sinc(C)  =  Pi  •  vx  +  P2  •  ^2  • 

Transforming  produces  the  filter  X[~p,0\. 
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